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Abstract 

A system with a metastable phase and a pseudo continuous set of 
the heterogeneous centers is considered. An analytical theory for ki- 
netics of the process of condensation in such a system is constructed. 
The free energy of formation of the critical embryo is assumed to be 
known in the capillary (macroscopic) approach as well as the solvata- 
tion energy of the embryo. The theory is based on the quasistationary 
approximation for the nucleation rate which has been justified analyt- 
ically. An effective iteration procedure is presented. The iterationa 
are calculated analytically. The approximate universal form of the 
spectrum is established. 
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1 Introduction 



The theory considered here completes the program of construction of the 
theory of the aerosol formation announced in [14|. It will be based on the 



capillary approximation of the height of the activation barrier. This approxi- 
mation lies in the base of the classical theory of the homogeneous nucleation. 
All necessary bibliographic remarks can be found in |TJ[ . 



Speaking about the activity of the heterogeneous centers we mean that 
the different activity initiates the different height of the activation barrier 
AF, i.e. the difference between the free energy of the critical embryo F c and 
the free energy of solvatation G. The set of the different activities of the 
heterogeneous centers can be so dense that we can regard it as the pseudo 
continuous one. In the case of the "solid nucleus of condensation with the 
weak interaction" the continuous size of the nucleus ensures the continuous 
set of AF. 

The ordinary external conditions for the phase transition have the smooth 
character in time. The external action on the system leads to formation 
of the droplets of the new phase. The process of condensation violates the 
thermodynamic parameters of the system. When the external action on these 
parameters has the smooth behavior in time we shall say that condensation 
occurs under the dynamic conditions. 

Nevertheless the theory of the heterogeneous condensation in the dynamic 
conditions has been constructed only for one type of the heterogeneous cen- 
ters. So, the task to construct the kinetic theory for the system with a 
continuous set of the heterogeneous centers is rather essential. It will be 
completely fulfilled here. 

We shall use the physical assumptions analogous to jnj which are neces- 
sary to construct the mathematical model: 

• the thermodynamic description of the critical embryo, 

• the random homogeneous space distribution of the heterogeneous cen- 
ters, 

• the free-molecular regime of the droplets growth, 

• the homogeneous external conditions for the temperature and for the 
pressure, 
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• rather a high activation barrier. 

As far as the most interesting characteristics of this process are the numbers 
of the heterogeneously formed droplets on the centers with a different activity 
we shall estimate the accuracy of the theory by the error of the obtained 
solutions for these values^. The unit volume is considered. All energy-like 
values are measured in the thermal units. 

The publication has the following structure: 

• In the first part the system of the equations of condensation is con- 
structed. 

• In the second part the explicit calculation of the iterations is presented. 

• In the third part the conception of the universal spectrum is developed. 

• In the last part some realistic spectrum of activities is considered and 
the principle of the self-consistency of the spectrum of the activities is 
presented. 

We shall define the activity of the heterogeneous center as some parameter 
w which is proportional to the height of the activation barrier 

AF(w) = AF \ w=0 -Xw (1) 

with some positive coefficient A. Note that the choice of the initial point 
w = is rather arbitrary now. We suppose AF(w) ^> 1 for all essential 
types of the heterogeneous centers^. 

The total number of the heterogeneous centers with the given activity w 
will be marked by r] to t{w). Naturally t] tot (w) is rather a smooth function of 
w. We shall suppose that r] tot in the essential region is near some constant 
value or it can be well approximated by the polinom which power isn't too 
high. This value is constant in time. 

The density of the molecules in the equilibrium vapor is marked by n^, 
the density of the molecules in the real vapor is marked by n. The power of 
the metastability will be characterized by the value of the supersaturation 



1 But not for the total number of the droplets. 

2 When this condition is violated the kinetics of the process can be described by the 
trivial modification of the theory presented below. 
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We shall define the super-critical embryos as the "droplets". Every droplet 
is described by the number of the molecules v , or by the linear size 

p = v 1 / 3 

Due to the free-molecular regime of the droplets growth we have 

dp _! 

Tt =CaT 

where a is the condensation coefficient and r is some characteristic time 
between the collisions in the saturated vapor obtained from the gas kinetic 
theory. 

Let us introduce the size z according to 

z = j\aT~ l dt' (2) 

Here £* is some characteristic moment of time which belongs to the period of 
the intensive formation of the droplets. The choice of £* is rather arbitrary. 
One can use for example the choice described in |J. Until the beginning of 
the coalescence 0,0 which isn't considered here equation (g) ensures the 
growth of z in time and can be inverted 

t(z)= ( Z ra- l -^- + U (3) 
Jo C{x) 

Hence, all values dependent on time become the values dependent on z and 
the relative size 

x = z — p 

can be introduced. During the whole evolution the droplet has one and the 
same value of the variable x. Considering t(x) as the moment when the 
droplet with the given x has been formed (as a droplet) we can consider all 
functions of time as the functions of x . Hence, we can see that the kinetic 
equation is reduced to the fact that every droplet keeps the constant value 
of x. To reconstruct the picture of the evolution one must establish the 
dependencies t(z) and ({x). 

The values at the moment £* will be marked by the lower index 
The positions of the region of the intensive formation of the droplets are 
essentially different. But one can introduce £* as the moment corresponding 
to the maximum af the supersaturation. 
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2 The system of the equations of condensa- 
tion 



We shall mark by the argument oo the total values of the magnitudes formed 
during the whole condensation process. 
Introduce the value of 

u _ ntot 1 

Sideal A 

where n to t is the total number of the molecules in the system. 

We must take into account the reduction of Qdeai to some value $ due to 
the consumption of the vapor molecules in the process of solvatation || [10 



PH . Moreover, according to jTTJ the ideal supersaturation can be changed 
by the external supersaturation. In any case we shall assume this value as 
the known one and mark it by $. In some rough approximation Qdeai ~ 

We shall mark by rf ot the total number of the heterogeneous centers of 
all types: 

rf ot = I dwr] tot (w) (4) 



where 7] to t(w) is the total number of the heterogeneous centers of the given 
type (the density of the total value of the heterogeneous centers). 
The following statements are valid in the further consideration: 

• (1) The main role in the vapor consumption during the evolution is 
played by the super-critical embryos, i.e. by the droplets. 

• (2) The quasistationary approximation for the nucleation rate is valid 
during the period of the essential formation of the droplets for those 
sorts of the heterogeneous centers which aren't completely exhausted 
in this process. 

The justification of the second statement^ uses the estimate for times if 
of the relaxation to the stationary state in the near-critical region which can 
be found in Jl], || (for the investigation of the heterogeneous barrier the 
consideration is the same one). 

3 The second statement isn't valid for those sorts of the heterogeneous centers which are 
going to be completely exhausted. But there the result is obvious. Note that the periods 
of the intensive formation of the droplets on the centers with the different activities don't 
coincide. 
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Let f s be the stationary value of the distribution of the sizes of the hetero- 
geneously formed droplets measured in the units of n^. It can be presented 
in the following form 

/• = fd(( x )> w )v(x, w ) ( 5 ) 

where rj(x, w) is the density on activities of the number of the heterogeneous 
centers which are free from the super-critical embryos and is given by the 
following formula M 

= W c + exp(-AF((,w))r 

where W + is the number of the molecules in the vapor which interact with the 
droplet in the unit of time, A e u is the width of the equilibrium distribution 

u=(u c +u e )/2 

A e u = J2 exp(- J P(i/) + G) 
and A c v is the halfwidth of the near-critical region 

2 l/2 

A r V = 



■ d 2 F \ 1 1/2 

, Q v 2 )U=U C | 



Index "c" marks the values for the critical embryo and "e" - the values for 
the equilibrium embryo. Certainly, A c z/ and A e z/ are some smooth functions 
of w and we shall neglect this dependence. 

We shall mark by n 00 g{w) the density of the distribution on activities 
of the total number of the molecules of the condensated substance in the 
heterogeneous droplets formed on the centers of the activity w. To simplify 
the formulas we shall use 

9{w) = r](w)/r] tot (w) 
We obtain for g iy 9i the following equations 

g(z,w)= [ (z - x) 3 f c (((x),w)ri(w)dx (7) 



6(z,w) = exp(-rioo / f< : (((x),w)dx) (8) 



oo 
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As far as we measure the accuracy of the theory in the terms of the error 
in the droplets number we define these values as the following ones: 



N(z,w)=r) tot (w)(l-6(z,w)) (9) 
The total number of the droplets is 

N tot = f r] tot (w)(l - 6(z,w))dw = f N(z,w)dw (10) 



For the majority of the types of the heterogeneous centers the following 
approximations of the nucleation rates are valid during the period of the 
essential formation of the droplets^] 

fc((( x )> w ) = w ) U=o exp(r — — - -)exp(wA) = f c *exp(r — — - -) exp 

(11) 

where 

dAF(Q 

i - -®* — ^ — lc=*.U=o 
fc * = fd®*' w ) U=o 

and AF is the height of the heterogeneous activation barrier^. The validity of 
these approximations can be justified for the heterogeneous embryos with the 
interaction between the center and the molecules of the condensated phase 
weaker or equal than the reciprocal to the space distance. So, we can imagine 
this as a hard sphere with a weak interaction on which the embryo is formed. 
The dependence of T on w is rather weak. So we can put 

i» = r u =0 (13) 

for any essential part of the spectrum of the activities. The applicability 
of the last approximation is based on the following qualitative model. It 

4 Note that the position of the regions of the intensive formation of the droplets on some 
sort essentially depends on the activity of the centers, the following approximation isn't 
valid for the regions of the intensive formation on the active centers. But for these centers 
the result is evident - all centers are going to become the centers of the droplets. 

5 The validity of these approximations may fail in the situations of the extremely long 
spcctrums (see later). But the square form of the behavior of the supersaturation near 
the maximum will be attained earlier. 
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is known that T is approximately $*(f c — v e )/(&* + 1). When the super- 
saturation is sligtly changed then the value of u e doesn't undergo some big 
variations and the variation of V occurs mainly due to the variation of v c . The 
behavior of the value v c resembles the analogous behavior in the homogeneous 
case. As a result we can justify not only the approximation concerning the 
dependence on the supersaturation, but also the dependence on the activity. 

A natural question about the essential part of the spectrum appears here. 
The process of condensation on the centers with some relatively high activ- 
ity occurs earlier than the supersaturation attains maximum. These centers 
form the droplets which can be included in the value of the external super- 
saturation in the manner of |~T| [[16]. For these centers the number of the 



formed droplets coincides approximately with the number of the heteroge- 
neous centers - all centers are now the centers of the droplets. 

The action of this part of the spectrum on the further evolution can be 
treated in the terms of the external supersaturation Q (see The non- 

trivial statement that the process of formation of the droplets on the centers 
with the intermediate activity occurs near the maximum of the supersatu- 
ration lies in the base of such a method of description. This fact can be 
analytically proved. 

The centers with some relatively low activity remain practically unex- 
hausted (when the spectrum f] tot has no singularities they can not play any 
important role in condensation). The singularities can be described in this 
region in a manner from |T| ||16|| . So only some centers with the intermediate 
activity are essential near the maximum of the supersaturation. In the scale 
of the activities this region corresponds to the variation of w by the value of 
the order of A -1 . 

Now we shall formulate the system of the balance equations for the func- 
tions g(x,w), 9(x,w) and 

G(x) = J dwg(x,w) 

where the integral is taken over the whole spectrum of the heterogeneous 
centers. In principle we can write that the region of integration goes from — oo 
up to oo but in fact this integration must be carried out in the region which 
covers the essential region of the spectrum^. In this region it is necessary to 

6 The extraction of this region is quite analogous to the case of the decay on the spectrum 
of activities. 
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have relatively smooth behavior of r] to t- This value must be approximately 
constant 

r] tot ps const 

or must be well approximated by the polinom on the activities: 

Vtot = Pn{w) 

We assume the total number of the heterogeneous centers to be constant 
in time. 

Using the conservation laws for the heterogeneous centers and for the 
molecules of the condensing substance we get for g, 9 the following equations 

(z — x) 3 exp(— r — — — -)9dxexp(w\) (14) 

-OO Mr ^ 

G(z) = J dwg(z,w) (15) 

at \ ( /*exp(Aw)ra 00 r z C~$*n. x na\ 

v(z,w)=exp( / exp(— r )dx) (16) 

Vtot J-oo 

$ = ( + G(z) (17) 

where /* = *n tot 

These equations form the closed system of the equations for condensation 
kinetics. This system will be the subject of our investigation. 

We shall consider this system of equations during the period when the 
centers with the intermediate activity are going to become the centers of the 
droplets. We shall call this period as the period of the intensive formation of 
the droplets (PIFD). 

We assume that the ordinary || linearization of the ideal supersaturation 
is valid during PIFDQ 

$(x) = + y cx ( 18 ) 

7 One can analytically show that the duration if PIFD is rather short. Namely 

$(z)-$. 1 

- r 
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with some positive parameter c. 

After the substitution of this linearization into the system of the conden- 
sation equations this system transfers to 

r z r 
g(z,w) = /* / (z — xf exp(er — — G(x))9(x, w)dx exp(w\) (19) 

J—oo $* 

G{z) = J dwg(z,w) (20) 

at \ t /* ex P( Aw ) n °Q f z i 1 m wa \ fot\ 
9(z,w) = exp( / exp(cx — — G(x))dx) (21) 

Vtot J-oo 

The spectrum of sizes can be found as the following one 

c - 

f(x,w) = /* exp(Aw) exp(— T — — — )9(x,w) (22) 
and when the linearization is taken into account: 

r 

f(x, w) = /* exp(Xw) exp(cx — —G(x))9(x, w) (23) 

3 Iteration procedure 

Such systems as the already obtained one can be solved by the iteration 
procedure. It can be constructed by the following way: For the initial ap- 
proximations we choose: 

g (z,w) = 9 {z,w) = l (24) 
The recurrent procedure is defined according to 

gi+i(z, w) = /* / (z — x) 3 exp(cx — —Gi(x))9i(x, w)dx exp(wA) (25) 

J—oo $* 

Gi{z) = J dwgi(z,w) (26) 



a ( \ < /*exp(Aw)n 0O f z ( V n i \\a \ tws 
9 i+ i(z, w) = exp( / exp(cx — —&i(x))dx) (27) 

Vtot J-oo <P* 
The chains of inequalities analogous to [fjj guarantees the convergence of 
the iterations and some estimates analogous to J7| can be established here. 
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The direct calculation of the iterations gives 



gi(z,w) = /, 



6 exp(c^) 



exp(Aw) 



(28) 




(29) 



(30) 



The last integral causes some problems with convergence. When rjtot = const 
it can not be calculated. At least there are two possibilities to overcome this 
difficulty. 

The first possibility is a more formal one. Certainly, in every system 
Vtot{w) = when w is greater than some maximum value w max . So, formally 
the integral exists and it is equal to some constant. This constant is generally 
unknown, but one can act in the manner like it was done with the supersat- 
uration in ||. In || the value of was not the matter of consideration in 
the iteration method, but was included in the set of the undefined parame- 
ters. After the final form of the size spectrum was established an equation on 
parameters was studied. The form of the spectrum can play a role of some 
ansatz with several parameters. The same procedure can be inserted here. 
The reason is that we can observe the separation of the functions containing 
w from the terms containing the other parameters of the condensation pro- 
cess. This fact remains valid also in the second iteration which is rather a 
good approximation for the final solution as it can be analytically proved on 
the base of the analysis of the iteration procedure. 

The second possibility seems to be more physical one because it allows 
to combine the current problem with the problem of the correct definition 
of the region of rather active and exhausted centers. Really we have not 
separated the region included into the external supersaturation from the 
region which isn't included into the external supersaturation. Note, that in 
the dynamic conditions the spectrum of sizes of the droplets formed on some 
fixed sort of the heterogeneous centers has well defined boundaries. The sizes 
(and the times of formation) of the droplets which are essential in the vapor 
consumption can be also well localized. Namely, the spectrum is essentially 
situated in the region 




(31) 
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and the region of the droplets essential for the vapor consumption during 
PIFD is essentially covered by the interval^ 

- ^ < z < (32) 

As a result one can state that for z ~ —a/c where a ~ 8 one can get the 
start of PIFD for the centers which are in the process of nucleation near the 
maximum of the supersaturation. 

Note that the moment can be taken as a moment of the maximum of 
the supersaturation. Really, when the maximum of intensity of formation (or 
the moment of formation of the half of the droplets) lies far from maximum 
of the supersaturation, then the result of nucleation of this sort of centers 
is rather obvious: the supersaturation is equal to the ideal (or external) 
supersaturation 

c = $ 

the total number of the droplets is equal to the total number of the hetero- 
geneous centers 

Ntot = Vtot 

and the moment of formation of the spectrum on the centers of some fixed 
sort which can be considered as a monodisperse one can be found from 

6{z,w) = 1/2 

which is equivalent to 

* exp(Aw)noo exp(cz) = cln2 

where the first two terms can be changed by some slightly other amplitude of 
the spectrum and the parameter c can be reconsidered as far as the derivative 
can vary. 

The simplest regularization by the cut-off is the following one. We must 
substitute the value of r\ tot by the value r\ init calculated as 

Vinit = VtotO(x = -a, w) 
8 See the analysis of the second figure from |lq| . 
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after the calculation of 6(—a,w) which can be calculated on the base of the 
ideal supersaturation, i.e. in the first iteration. As far as 

of \ , , (x ^exp(-ca) 
9i(-a,w) = exp(-/ c ^n^ exp(Aw) ) 

the value of f\ must be changed by 

J* = ft *Vtot{w) exp(-/ c exp(Aiuj ) 



4 Calculation of the iterations 

Let us calculate the iterations. Then the expression for G will be the following 

6 exp(cz) 



c 4 



f°° ( 
G\ = ft * / ^tot(w) exp(Aw) exp(-Aexp(Aw))(iM;- 

where 

A = k , n j^l 

When rj tot = const the integral can be taken which leads to 

1 6exp(cz) 
Cti - /c *^^4 - 4 

Note that the integral can be taken in the final limits which allows to take into 
account only the finite region of activities as far as an approximation of the 
total number of heterogeneous centers by the composition of the Heavisaid's 
functions. 

When r/tot is approximated by the polinom it is easy to note that each 
monom in the integral leads to some elementary functions with some uni- 
versal constants. This fact can be seen from the simple translation of the 
variable in _ 

w m exp(Xw) exp(—Aexp(\w))dw 

-00 

to the variable 

Xw + \n(A) 

Here the integration must be fulfilled in the infinite limits. These integrals 
can be expressed through \l/-function and its derivatives when the argument 
is put to 1. 



13 



The second approximation for 9 leads to the following result: 
9 2 (z,w) = exp(— f(_ ^n OD exp(Xw) f exp(cx — B exp(cx))dx) 

J —oo 

where 

T 6r] tot 1 

B = irk 



$*■"* c 4 XA 

Note that the integral can be taken here also in the finite limits. After the 
integration we get 

9 2 (z,w) = exp(-/ c „n 00 exp(Aw)^-(l - exp(-S exp(cz)))) 

CD 

Note that r\{w) = r] in i t (w)9(z,w). 

Particularly, for the final values one can get 

6> 2 (oo, w) = exp(-/ c ^noo exp(Aw)- 1 -) (33) 

CJD 

This value of B is obtained in the approximation rj = const. For the fur- 
ther consideration only the analytical structure of B is important. Namely, in 
any case (and also in the case rj tot {w) = P n (w)) we have for G the expression 
with the following analytical structure 

Gi ~ const(z, w) exp(cz) 

So, B is some constant and this fact ensures the possibility of the further 
calculation of the iterations. 

The next step is the calculation of the total number of the droplets ap- 
peared in the considered period which can be done by the simple integration 

/oo 
dwN tot (w) 
-oo 

Note that instead of 

/oo 
dwr kot (w)(l-9 2 (< X >,w)) 
-oo 

which can not be integrated we must take 

/oo 
dwrj init (w)(l-9 2 (oo,w)) 
-oo 
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which gives the integral with no problems of convergence. To calculate the 
last integral note that 

r-a y 

02(—a,w) = exp(— fa * exp(\w)n OCJ / exp(ca; — — G)dx) 

J —CO 

Consider now the case r] tot = const. After the evident renormalization of 
the variables of integration the integral can be reduced to Jf^ exp(— exp(x)) (1- 
exp(— H exp(x)))dx where H is some constant. This value can be approxi- 
mately calculated by the following procedure. Decompose the internal expo- 
nents and get 



f 

J —i 



exp(— exp(x))(l — exp(— H exp(x)))dx 



/ 



i i 

ft(! 7-Tm)dx = (34) 



_ m v^oo cxp(ix) ^poo 



v oo 

^i=0 (j+1)! 



n y^oo Ml y_L 

Note, that it is necessary to left the serial in the denominator. 

The necessary accuracy will be ensured by the first three terms in the 
decompositions of 

Y^oo y^2 WP+i 

[°° ^j=o W d „ f°° (i+lF d 

Then an integral from the rational function can be simply calculated (also 
in the finite limits). 
Now the integrals 

V oo y l H l+l V n-1 y l W+ l 

T ~ [°°Wv i= ° (m)! dv~ (°°Wv ° (m)! dv 

will be calculated. These integrals appear when the polynomial approxima- 
tion for rj to t is accepted. 
Consider the function 
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as the function of a complex variable. Integrate this function along the closed 
path Q constructed in the following manner: the big circumference with an 
infinite radius; two straight lines ]y + iO, oo + i0[, ]y — iO, oo — i0[; the little 
circumference with a zero radius. The integrals along the circumferences are 
going to the zero. The integrals along the straight lines give: 



/ K+' y - (In, 4- tofl ^.g ^ 



which can be reduced to Ij and J2k<j e klk with some known coefficients. On 
the other hand this integral can be reduced to 



E 



1 y z H i+1 
i=0 (i+1)! 



E res(\n j+1 y — * ^ ) 



inside Q. As the result we have the recurrent procedure which allows us to 
calculate all Ij. 

Now let us see how the corrections due to the cut-off can be introduced 
in the value of g(z,w). In the first approximation we have: 

gi(z, w) = exp(Aw)(exp(cz) — exp(— ca)) 
The value of G\ in the approximation r) tot = const is given by 

Gi = k *~^^x^ ex ^ cz "> ~ ex P(~ m )) 

Note that it can be calculated for the finite band of the spectrum of the 
activities. The calculation under the polinomial approximation for rjtot can 
be done in the same manner. 

Instead of the pre-exponential factor one can put some other function. 
The separation of this expression into two factors depended on w and on z 
is obvious. 

For 6 2 we get: 

9 2 (z : w) = exp(— fa *nooexp(Aw) / exp(ca; — -B(exp(cx) — exp(— ca)))dx) 

J —oo 
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and after the simple calculation: 



6 2 {z,w) = exp(— <t n 00 exp(Aw)exp(i?exp(— ca))— — (1 — exp(— £? exp(cz)))) 
The final value is the following: 

8 2 (z, oo) = exp(— /<- i(; ?T, 00 exp(Au7) exp(5 exp(-ca))- 1 -) 

CD 



It differs from ( p3|) only by the remormalisation of the amplitude thanks to 
exp(£> exp(— ca)). 

The value of N tot is calculated in the same manner. 

Note that one can change — oo to —a and repeat all calculations in the 
same manner. 

The analytical estimates show that the second iteration gives rather a 
good approximation for the process. The reason is quite similar to || . Really, 
one can simply integrate all estimates from || to justify the validity of the 
obtained expressions. 

Another important remark must be given. One can see that all obtained 
values have the separation of the expressions (or arguments of exponents) in 
two terms: the first one depended on w and the second one depended on z. 
This reduction can not be observed in all high iterations. It is responsible 
for the absence of the cross influence of exhaustion of the heterogeneous 
centers with the different activity in the first iterations. The similar property 



is absent in the situation of decay [14|, [15|. Here such a cross influence 



can be neglected because the equation on the parameters of condensation 
balances the time of formation with the characteristic intensity of formation 
and leads to the fact that the characteristic width of the size spectrum is 
always approximately equal to -. The probability to form the droplet on the 
center is always determined only by the supersaturation and by the value of 
activity (and not by the total number of centers of the given sort). 



5 Universal solution 

Now the problem of the construction of the universal solution will be con- 
sidered. The idea of the universal solution || lies in the fact that after for- 
mation of the spectrum of the droplets the further evolution of the process 
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depends only from the first three (and zero) momentums of the distribution 
function. So if one can choose some variables in which the solution (the dis- 
tribution function) is the universal function (undepended on the parameters 
of the problem) then the expressions for the momentums have rather a sim- 
ple analytical structure combining some parameters of the problem with the 
universal constants. 

For the process of the homogeneous condensation the universal solution 
was obtained in 0. For the process of the heterogeneous condensation the 
universal solution is absent, but some pseudo universal solution can be con- 
sidered as the base for the further iterations |7]]. 

We shall rewrite the system of the condensation equations in the terms 
of C and $ 



+ -^cz = ( + G 

r 



G 





* 



— oo 



Introduce the function 




Then after the substitutions 



Xw 



w 



X 




G 



r 



9 



9 



one can get 



-z = 6 + G 
V 



G 
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d( w ) = ~w?T — I / {z — x) 3 exp(5)9(w,x)dxexp(w) 
AQ^p J— oo 

6(w) = exp(—fc exp(w) / exp(S)dx) 

V J-oo 



Let us choose p as 



V 



1 



and now there are no parameters in the expression for g. If the moment of 
£* is chosen as the moment when the supersaturation attains the maximum, 
then 

dw3 I (z — x) 2 exp(5)8(w, x)dxexp(w) 



p 

Then 

nop _ fj'n^l^l" 

k * p m*r 1/A 

An ordinary and natural condition to establish the zero point of activity 
can be written as the following one: 

9(w = 0, z = oo) = - 

which gives 

, ^oo _ In 2 

and there are no parameters in the last equation. 

As the result there are no parameters in the system of the condensation 
equations and the solution has the universal form. It can be analytically 
proved that it is a unique solution. All consequences now coincides with 
analogous conclusions from [[| 

Note that there is a property of a very smooth dependence of the number 
of droplets on the /*. So we can get the equation on the parameters of the 
process in some rough approximation (see the iteration procedure) and then 
use the universal law. 

The last model has the evident disadvantage. The equation for G demon- 
strates no convergence. Really, we can substitute the expression for g into 
the expression for G and fulfil the integration over w. Then we come to 

G = A 00 (z- xf- U , w , dx 
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where 

The evident necessary asympote G — > when z — > — oo leads to 5(z) — > ^ 
when 2: — > — 00. Then the exponents can be cancelled and the subintegral 
expression has the asymptote (z — x) 3 which doesn't ensure the convergence. 
Then we need some regularizations. The most evident one is to notice that 
after the process of any appearing of the droplets the constant value of rj 
transforms into 

Vnew = Void exp(-Q exp(iu)) 

where Q is some constant. 

We shall use this value as the initial one and after the analogous trans- 
formations we can come to 

G = J^-( Z {z- xf _ eXp( ^ } dx 

J — 00 



A$*p 4 7-oo V 1 Q + ^ f_ oo zexp(6(x'))dx' 

The coefficient can be made equal to some constant (may be 1) due to 
the choice of the scale of activities. This choice states some concrete value 
of Q (which has the same dependence on the choice as ^*^°° has). The limit 
Q — > corresponds to the already observed situation. Now we have more 
general equation taking into account the power of the previous deformation 
of the activity spectrum. Evidently the solution depends on the parameter Q 
and isn't universal. We see the further inclusion of the distribution into some 
more general set of the solutions for the various powers of the deformations 
of the spectrums by the previous nucleation. 

Due to the problems of convergence for Q = we must describe the form 
of the spectrum explicitly. We are going to show that the size spectrum has 
the universal character. 

When Q is going to zero the region where Q is unessential starts at z\ 
corresponding to 

r*i 

/ exp(S(x))dx ~ (p/c) exp(czifp) ~ Q 

J —00 

and becomes larger. The asymptote (z — x) 3 of the subintegral expression 
is attained from z% till some z with no respect to Q. So, the asymptotic 
region becomes larger. Certainly, we can integrate the asymptote explicitly 
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and come to (z — zi) 4 /4. Only due to the linear size this asymptote provides 
the main quantity of the substance in the droplets. As a result we have the 
following expression for the behavior of 



f c ~ exp(^ - a(z - ziY) 



where the scale of z is chosen to put c/p = 1 and a is some constant. The 
value of a can be determined from the requirement that max * (the max- 
imum of the supersaturation) is attained at z = 0. Then 

a = -zf 3 /A 

The expression for * can be presented in the following form 

* ~ exp(— az l ) exp(— (z/A) 2 ) exp(— 4o;2; 3 |2;;|) exp(— az 4 ) 
where A is the characteristic halfwidth 



One can see that at the characteristic scale z ~ A the arguments of the 
third and the forth exponents have the characteristic values z l and ir- 
respectively. Thus, these terms can be neglected and the form of * is the 
gaussian one: 



that the universal gaussian form of exp(5) ensures the universal expression 



This universal form will be attained earlier then Q = (where the property 
of convergence disappears) and has the true physical sense. 

Figure 1 shows the forms of the size spectrums in some different situations 
presented in the normalized coordinates. The curves "a", "b", "c", "d" are 
the solutions of the following equation 



A = (2|z z |/3) 



1/2 




for 9 



9 = exp(— exp(w)—== / exp(— x 2 )dx) 




z = Ai[ (z- xf 

J — oo 




dx 



This rescaling excludes the unphysical parameter z/. 
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for the different values of Q ( "a" : Q = 10 "b" : Q = 1 "c" : Q = 0.001 
"d" : Q = 0.0001 ). The curve "a" practically coincides with "b", the curve 
"c" practically coincides with "d" for all x. The curve "e" is the spectrum of 
the gaussian type exp(— x 2 ). The curve "f is the universal form of the first 
iteration exp(a; — exp(x)) after the renormalization in order to have the same 
position of the maximum of the curve and to conserve the integral over the 
spectrum. This curve practically coincides with "a", "b" but only for the 
negative values of the argument. 

One can prove that for all possible values of Q the spectrum lies between 
the first iteration and the Gaussian spectrum. It can be seen that the differ- 
ence between the first iteration and the Gaussian curve is rather small (less 
than 0.2) and lies in the frames of the accuracy of the modern experiment^]. 
Note that the value of Ai must be chosen to ensure the position of the max- 
imum of the supersaturation near max 5 = 0. For Q — > the quantity of 
the substance in the droplets isn't small due to the known power asymptote 
and one have to remove the choice of t* from the previous condition that 
the max 5 is attained at z = to the new condition max 5 = 0. In the 
concrete calculations one needn't to fulfil this condition precisely but only in 
a very approximate way as far as the form of the spectrum doesn't change 
essentially for the different choice of the max 5 near zero. This property can 
be shown analytically and is illustrated by Figure 2. 

Figure 2 demonstrates the weak dependence of the form of spectrum on 
the choice of the base point for the linearizations, i.e. the choice of t*. This 
property is important for the justification of the conception of the universal 
spectrum and can be proved analytically. This property is also necessary for 
the simplification of the numerical calculations. For Q = 10 the two values 
of Ai have been chosen. This corresponds to the different manners of the 
choice of or the choice of max 5. Two curves for 5 for different A\ are 
shown in the part "a" of this illustration. The upper curve corresponds to 
A\ = 0.5, the lower curve - to A\ = 1. In the part "b" the size spectrums 
for these situations are drawn in the normalized coordinates. They coincide 
and one can see only one curve drawn in the part " b" . 

Note that the limit Q — > is analogous to the situation of the "wide spec- 
trum" considered in Here we have some additional difficulties. These 

10 This ensures the weak sensitivity to the choice of the parameter a. It is necessary to 
the self-consistency of the presented theory. 
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difficulties appear from the asymptote (z — x) 3 for the subintegral expression. 
As a result the relative quantity of the substance in the droplets is greater 
than T" 1 , the spectrum is going to be formed during the period with duration 
essentially longer than c _1 r/£, the quasistationary approximation isn't valid 
during all the period of the droplets formation. So, one has to reconsider the 
approximations which form the base of the concrete constructions. 

The simplification goes from the asymptote (z — x) 3 . Really, this asymp- 
tote doesn't depend on the concrete form of the free energy and the deriva- 
tives of the free energy. So, the quasistationarity and exponential approxima- 
tion (|TT| ) aren't so essential (they are going to fail namely at this asymptote). 
The big quantity of the substance in the droplets require to reconsider the 
condition of the choice of t* as it is done above. The essential duration of 
the period of the droplets formation may lead to some new approximations 
and, thus, change the form of the asymptote (z — x) 3 ■ But the consideration 
of the process remains quite analogous to the already fulfilled one. One of 
such examples is given in the next section. 

The opposite situation corresponds to the validity of the iteration pro- 
cedure with the ideal supersaturation taken as the external one[5 Note 
that the final results for 9 and for the spectrum are obtained (in the sec- 
ond approximation) on the base of the supersaturation calculated without 
the appropriate account of the exhaustion of the centers, i.e. in the pseu- 
dohomogeneous situation. The same situation occurs when we study the 
condensation on the separate sort of the heterogeneous centers |§, @. This 
analogy leads to the three important consequences. To establish the first 
one note that the "pseudohomogeneous" base for the final results evidently 
has some universiality. In the universiality can be broken by the change 
of the regime of the consumption of the vapor by the droplets. The result 
of H is more approximate, but it is based on the form of the first iteration 
which can not be depended on the regime of the vapor consumption. So, the 
universiality here is more general. The analogy of the forms of the "pseudo- 
homogeneous" spectrums allows to give the definition of the relative activity 
also to the centers in the the process of condensation on the centers of the 
separate sort. The point w = corresponds to Ofi na i = 1/2. The value of w 

11 Certainly, with the corresponding choice of the renormalization oiOmu and the choice 
of the "external" supersaturation the iteration procedure remains valid in all cases with 
the linear "external" supersaturation. 
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can be reconstructed from 9fi na i as 



= ln(^4) 
l ln(l/2) J 

The invariant ln#/ exp(u>) can be also observed for an arbitrary z. Hence, the 
condensation occurs in a hierarchical manne iQ. Then one can consider the 
distribution f(x, w) over the two variables x and w obtained in the process of 
condensation on the heterogeneous centers with the spectrum of activities as 
the formal generalization of the distribution in the process of condensation on 
the heterogeneous centers of the separate sort. Now the distribution has the 
universal functions as the base. This base doesn't depend on the process of 
the exhaustion of the heterogeneous centers. The dependence on w is rather 
explicit 

rx \ n 2 
f(x, w) ~ exp(5(x)) exp(— w / exp(5(x'))dx' - 



J^ 00 exp(5(x'))dx' 

The factor exp(<5) is the universal base for this distribution. That's why we 
prefer to present it in Figure 1 and to call it simply the spectrum. The 
last expression is precisely valid for 5 calculated on the base of the precise 
solution. As for approximations for S one can take 5 from the pseudoho- 
mogeneous situation in two manners (from |5| and from 0). The first one 
corresponds to the first iteration and doesn't depend on the regime of the 
vapor consumption. The second one is more precise but depends on the 
regime of the vapor consumption. Note that the last expression contains 
an unessential approximation because it is based on the pseudohomogeneous 
situation. The negligible character of correction can be seen simple from 
the result of the iteration procedure in the situation with one sort of the 
heterogeneous centers. Hence, one can state the natural inclusion of the spe- 
cific condensation on a sort of the heterogeneous centers into the practically 
universal distribution appeared from the process of the condensation on the 
heterogeneous centers with the spectrum of activities. 

The second consequence is connected with the fact that now we know the 
ansatz necessary for the realization of the first opportunity (alternative to 

12 Certainly, the shift and the zero roint w = depend on the parameters of the process 
including r] tot . 
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the regularization of 8). One can use the form of 6 obtained in the pseudoho- 
mogeneous situation in the first iteration or from the universal homogeneous 
precise solution. 

The third consequence can be noticed in frames of the modified method 
of the steepens descent ||16|| . Really, the form of the external supersaturation 
can be essentially nonlinear^. Then we have to use the modified method of 
the steepens descent. The generalization is evident. Due to the negligible 
effect of the account of the exhaustion of the heterogeneous centers in all 
steps except the final formulas which can be seen from the iteration method 
(in the nonlinear case all analogous estimates can be given) we can simply 
integrate (summarize) the amplitudes of the distribution^ over the activity, 
then get the result for the psudohomogeneous situation and then apply the 
final formulas for the number of the heterogeneous centers on the base of the 
supersaturation. 

One can analytically show for Q > 1 that the differences in the forms of 
the spectrums between the real solution and the first iteration, between the 
real solution and the pseudohomogeneous universal precise solution decrease 
when Q increases. Also one can prove that for Q > 1 the difference between 
the external supersaturation and the real supersaturation at some arbitrary 
moment of time decreases^ 5 ! when Q increases. So, the results of Figure 1 show 
that the required limit corresponding to the practically universal spectrum 
is attained already when Q > 1. 

One can analytically show that the form of the spectrum lies between 
the gaussian iteration and the first iteration (or between the gaussian form 
and the universal precise pseudohomogeneous solution). We can see from 
Figure 1 that the difference in the form of the spectrum is rather small for all 
situations. This can lead us toward the general approximate universiality. So, 
the analogy with the ideas of the universiality in the previous considerations 
is now stated. 

Certainly, in the case of the condensation on the separate sort of the 
heterogeneous centers the spectrum (as 8exp(5)) isn't universal. But the 

13 In a very specific global form of the activity spectrum r] to t- One can prove that in any 
case one can restrict the expression for the external supersaturation during the period of 
the intensive formation of the droplets on the centers with the intermediate activity by 
the constant, linear and square terms. 

14 The characteristic values of the distribution at t*. 

15 And also at t*. 
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difference between the forms of the spectrum isn't so important. One can 
show analytically that the formj^] of the spectrum changes continiously from 
the case of the relatively small number of the heterogeneous centers up to 
the pseudohomogeneous case and lies between these two limit cases. Really, 
when the relative number of the heterogeneous centers is small the first it- 
eration exp(x — exp(x)) gives the precise expression for the spectrum. In 
the pseudohomogeneous case the same first iteration is the base for the re- 
sult which lies prctically near the precise solution (this is the base point for 
the effectiveness of the iteration method |§). So, the form of the spectrum 
doesn't essentially variate. It can be also treated as as some universal form. 
Hence, this case also allows the universal description. Namely this is the base 
for the applicability of the iteration procedures from ||, (J7|. 



6 Advantages of the model 

The results of the iteration procedure can give not only the description of 
the process of the heterogeneous condensation but also the information about 
the form of the activity spectrum. From the first point of view this question 
appears to be an external one to the process of condensation. 



Experimental results [ 1 3 1 show that the spectrums of the activities are 



rather smooth and have the form 

(p(v) ~ v- (1+s) 

where s is some small positive parameter. Here v has the sense of the activity 
introduced in a slightly another manner. Namely, the free energy can be 
approximately written as the function of v in the leading term as 

F = — ln(C + l)p 3 + constp 2 + constv In v + const 

which leads to the approximate applicability of the linearization of AF = 
F c — G = F(p c ) — F(p e ) as function of w together with the identification of 
v as w — w with some characteristic value of the parameter w . Then r\ tot 
can be presented as 

r]tot{w) = rj tot (w ){w - w o y {1+s) (35) 



16 After renormalization and as the function of a rescaled variable. 
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The positive value of s ensures the convergence of the total number of the 
heterogeneous centers^ 



V = / Vtot(w)dw 

J w 

Define the class of the "long tail spectrums" (LTS) as the spectrums with 
extremely long tails in the active region. This spectrum shows an example 
of such a spectrum. Here we shall develop the theory for condensation on 
such spectrums. 

To present the most simple variant we shall assume that for an ap- 
proximation of the known form is adopted 

r 

/c. = /c*» ex P(^(C* - C**)) (36) 

with some constant value of the known parameter T. Here can be treated 
as a value corresponding to the moment when some sort of the heterogeneous 
centers is exhausted and £** is the base for approximations. Require that 

9(z = 0,w) = 1/2 

which leads to 

7— ((* - C**) + Aw = In 



C** " cln2 
The last approximation is valid when the relative variation of ( is small. In 
such a region we can put c in the r.h.s. of the last equation to some constant 
value and get for the size spectrum 

I f( Z ) 1 = 1 Vtot(w ly.- const-cz ) I - 

x A 

where the dependence lnc is assumed to be unessential. Then one can see 
that in the case of (^) the convergence of the integral for G can be attained 
only if s > 3 (note that c — > const ~ ^ when w — > 00). This result 
shows that LTS can not be effectively spread. More rigorously speaking, 
the contradiction can be overcome by noticing that this effect is only due to 
approximation (p6|). Certainly, G is restricted by the value r/ tot ^ ax where 



1 The infinite tail must be cut off at the centers which must be already exhausted in 
the previous processes. 
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z max is the coordinate of the droplets which are imaginary formed when the 
supersaturation attains some slightly positive value. Nevertheless in such a 
situation all global features of the free energy appear which doesn't allow 
to get any effective method to all types of the heterogeneous centers. The 
leading idea will be that some of heterogeneous centers had been exhausted 
in the previous processes of condensation occurred earlier. Then the balance 
equation will be the following 



- = c- + SJ Jz-x) f(x)dx 



where d is some boundary parameter of size spectrum initiated by boundary 
of activity spectrum. This equation can be solved. In the case of r] tot ~ const 
we have 

(** . o [* I \ 2 / \ J COnst f O r 7\ 

■ = c— + o / [z — x) c{x)dx — - — (37) 
dx T J-d A 

Note that here c^fr is the derivative of the real (not the ideal) supersaturation 
on x. Having differentiated this equation three times^ we get the ordinary 
linear differential equation with some constant coefficients, which has the 
known solution. 

Now a new principle of the self-consistency of the activity spectrum will 
be elaborated. We have noticed that active centers have been exhausted in 
the previous processes of condensation. But the description of the previous 
processes of condensation is quite analogous to the description of the given 
process. As the result we have for spectrum of the heterogeneous centers 

9 = exp(— const exp(Au>) / exp(— T )dx) 

J — OG 

with the analytical structure 

9(z oo) = exp(— const exp( \w)) (38) 

But this structure is already known - it is the structure of the final and the 
start results for 9 in the already investigated process. So, our start form was 
absolutely correct^. Also it is very important that the structures before and 

18 Now we are going to get the equation for the unknown function c(x) so we have to 
keep the derivatives of c(x) 

19 One can repeat all calculations with the new value of parameters initiated by the 
previous process of consumption of the heterogeneous centers. 
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after the process coincide in their analytical form. This statement will be 
called as the principle of the self-consistency of the spectrum of the activities. 

Hence, our initial conditions satisfies this principle. The process of con- 
densation conserves the analytical form of the spectrum of the activities and 
the result of the process can be regarded as some shift in the scale of activities. 
Only some parameters are changed. We needn't to repeat the calculations 
as far as they are absolutely the same ones. 

However, in fl37|) when d > 1 we needn't to know the details of the cut-off 
of the spectrum and can use the cut-off by the Heavisaid's function. 

The last point of discussion is the possibility of the linearization of the 
external supersaturation. One can analytically prove the following statement: 

• In the description of condensation on the centers of the intermediate 
activity the two manners of description cover all possible situations^. 
These manners of description are the following: 

1) . The supersaturation is absolutely determined by the consumption 
of the vapor by the droplets formed on the centers with high activity. 

2) . The linearization of the external supersaturation is possible dur- 
ing the period of formation of the droplets on the centers with the 
intermediate activity. 

The first manner of description is obviously trivial. The second is com- 
pletely described here. Note that in the case when the spectrum has some 
pseudo singularities the last statement isn't valid and the special description 
is necessary. This description can be attained by the combination of the 



methods presented in [16| with this theory 



7 Concluding remarks 

Now the theory of aerosol formation is completed. As a base for some concrete 
results the classical theory of nucleation was chosen. Meanwhile the validity 
of this theory remains the matter of discussion^. Note that this theory is 

20 The natural requirement for the activity spectrum is that the activity spectrum must 
be a smooth function with the cut-off (or without the cut-off) like (Sq). 

21 Note that formally we can put the discrepancy between the predictions of the classical 
theory and the real rate of nucleation into the values of the microscopic corrections. 
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necessary only to calculate the amplitude of the distribution. In the situation 
of the dynamic conditions the manner of the choice of £* leads to the algebraic 
equation of the following form: 

/* = Smooth known function^*; explicit external parameters) 

So, the microscopic corrections to the free energy aren't important for the 
process of condensation under the dynamic conditions. In the situation of 
the metastable phase decay the amplitude value of the spectrum is rather 
artificial. It is given explicitly by the initial supersaturation and one can not 
exclude the sharp dependence on the microscopic corrections to the free en- 
ergy. Here any other concrete formula can be used instead of the expression 
from the classical theory for the nucleation rate in the stationary approxima- 
tion. Note that the form of the spectrum doesn't essentially depends on the 
value of the amplitude. That's why we don't pay any serious attention to 
the solution of the algebraic equations on the parametrs of the process and 
to the calculation of the concrete numerical results. We have concentrated 
our attention on the universal depencies which form the base for the common 
knowledge in the field of the first order phase transitions. Namely the univer- 
siality is the main result of our investigation. Note that the universiality also 
takes place in the process of the condensation under the dynamic conditions. 

The theory presented in these four publications can be easily reformu- 
lated when we assume that the rate of appearance of the critical embryos is 
proportional to 

F ie xp(F 2 ) 

where F\ and F2 are the smooth functions of the supersaturation and F2 has 
the big absolute values. Evidently all reasonable recipes for the rate of nu- 
cleation satisfyes the last ansatz. We kept the classical stationary flow only 
due to the tradition. The range of applicability of the presented theory is 
wider than the case of condensation. The possibility of the linearization of 
the ideal power of metastability during the period of the intensive formation 
of the super-critical embryos seems to be quite natural. The sharp increasing 
of the intensity of the metastable phase consumption by the separate super- 
critical embryo is also a natural requirement. This theory is based only upon 
these assumptions. All positions used in the construction of the mathemati- 
cal model (exept the absence of coalescence) can be missedP] by the obvious 

22 The thermal effects can be taken into account in the manner analogous to M. 
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generalization conserving the mathematical structure. The requirement of 
the homogeneous character of the external action and of the nucleus distri- 
bution can be attained by the consideration of an arbitrary hydrodynamic 
element. As the result one can conclude that the general theory for the ki- 
netics of formation of the dispersed embryos of a new phase on the dispersed 
unpurities in a metastable phase is constructed. 
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Figure 1 



Form of spectrum under the different powers of deformation. 
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Figure 2, part a 



1 

Figure 2, part b 



Weak dependence of the spectrum on the base of approximations. 
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